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Abstract 


Introduction. The development of the polar areas of the World Ocean and the need to solve various problems 
associated with a large number of freezing inland water bodies issue new challenges for science. These challenges 
include the problem of studying the behavior of ice cover when exposed to various types of loads. Of great interest is 
the consideration of problems about the action of a moving load on the ice cover. A moving load simulates the effect 
of moving vehicles on ice. However, in papers devoted to the above problems, cases of load movement along a straight- 
line trajectory are considered. The objective of this research is to develop a method for studying the behavior of ice 
cover under the action of a load moving arbitrarily. 

Materials and Methods. The article proposes a method for solving the problem of the action of a force moving along 
an arbitrary trajectory on the ice cover of a reservoir of finite depth. The problem amounts to solving a system of two 
differential equations. The first of them models the behavior of the ice cover, and it is the equation of vibrations of a 
viscoelastic plate. The second equation simulates the behavior of fluid in a state of potential flow, and it is Laplace's 
equation. To solve the system of differential equations, integral transformations in time, space and variables were used. 
The resulting solution was expressed through an iterated integral, which was calculated using numerical methods. 
Results. The development and implementation of the method resulted in solving the problem of the movement of a 
concentrated force along an ice cover according to an arbitrary law. At the same time, studies were carried out on the 
behavior of displacements and stresses in the ice cover depending on the speed and acceleration of the movement of 
the vertical load, on the depth of the reservoir, and on the viscoelastic properties of ice. In addition, the distribution of 
the velocity vector of fluid particles along the depth of the reservoir was calculated. 

Discussion and Conclusion. The proposed method is very effective for solving problems of moving loads acting on 
the ice cover of a reservoir of finite depth. It provides solving problems about the action of a load moving along an ice 
cover along a complex trajectory. The results obtained can be used to calculate the stress and displacement of the ice 
cover during the laying of ice roads or the construction of airfields on the ice. 
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AHHOTauna 

Beedenue. Ocsoenue nouApHbIx pationos Muposoro okeaHa, HEOOXOAMMOCTb pellieHHA pa3ssIM4HBIX 3aja4, CBA3AHHBIX 
C HaMyveM OONbWIOTO YMCA 3aMep3aIOWHX BHYTPCHHHX BOJOCMOB, CTAaBAT Tepe HayKOl HOBbIe mpoOsempl. K ux 
YHCIIy OTHOCHTcA MIpoOIeMa H3YHCHHA NOBeACHHA JeAAHOTO MOKpoBa MOA BO3AeHCTBHeM Ha Hero pa3IM4HOrO Bua 
Harpy30K. bombo HATepec NpeAcTaBsIAeT paccMOTpeHHe 3afay O elCTBHM Ha JIeqAHOM MOKPOB NOABWKHOHM Harpy3Ku. 
TloqBpywxKHad Harpy3ka MOjjeMpyeT AelicTBue Ha JIeq ABIMDKYIMXCA TpaHCOpTHbIx cpeycTB. OyHako B paodoTax, 
TIOCBAL[CHHBIX BbIMIeyKa3aHHbIM 3alayaM, paCCMaTpHBaIOTCA CilyyaH BMWKeHHA Harpy3KH MO UpaAMOMMHeHOH 
Tpaextopuu. Llenbro anno paooTsl ABIAeTCA pa3spaOoTKa MeTOJa HCCIeAOBAaHHA MOBeACHHA JIeAAHOTO MOKpoBa Toy 
TeHCTBHeM Harpy3KH, MepeMellaroulelica IPpOM3BOJIbHEIM OOpa30M. 

Mamepuaavi u memoool. B cratbe ipeqoxKeH MeTOX pelleHHA 3aa4uH O AevCTBHM Ha JIeqAHOM MOKpOB BOoeMa 
KOHeYHON PIYOHHEI ABMOKYWeHCA 10 IPOM3BOJIbHOM TpaeKTOpHH CuIbI. 3ayaya CBOAMTCA K PeIeHHIO CHCTeMBI JByX 
TuddepenunanbHEix ypaBHennn. Ileppoe 13 HUX MOJeMpyeT MOBeeHHe eqAHOTO TOKPOBa H ABJIACTCA ypaBHeHHeM 
KoseOaHHi BA3KOYUpyrou WiacTHHbI. Bropoe — MoyesMpyeT NOBeeHHe 2%KHAKOCTH, HaxOAulelica B COCTOAHHU 
MOTCHUMaIbHOrO TeYeHHA, MH ABIAeTCA ypaBHeHuem Jlanmaca. JIna pewleHua cHcTeMbI AUdpdepenuMabHBIx 
ypaBHeHHM TIPHMeHAIMCb WMHTerpasibHble MpeoOpa30BaHHA MO BpeMeHHOM HW MpOCTpaHCTBeHHbIM TIepeCMeCHHBIM. 
IlomyueHHoe B pe3yIbTaTe pelieHHe BbIpaxkaocb Yepe3 MOBTOPHBIM UHTerpall, AIA BbIYMCIeHHA KOTOPOTO 
IIPHMCHAJIMCb YACICHHBIC METOBI. 

Pe3yismamoet ucciedoeanua. B pe3yibTate peaM3alMn WpesoKeHHOrO MeTOJa HoMyYeHO pelleHve 3aqaun oO 
TBWKCHHH COCpexOTOYCHHOM CHJIbI 10 JIeqAHOMy MHOKpoBy MO MpOu3sBOIbHOMy 3aKoHy. IIpu 3TOM Mpou3BeeHbI 
McCCIeqOBaHHA XapakTepa MOBeCHHA MepeMeICHHH M HallpsKeHHH B JIe{HOM MOKPOBE B 3ABHCHMOCTH OT CKOpOCTH H 
YCKOpeHHA JBMOKCHHA BeEPTHKAIbHOM Harpy3KH, TyYOMHbI BOOeCMa UM BASKOYUPyrux CBOMCTB JIbya. Kpome toro, 
paccunTaHo paciipeyjeeHve BeKTOpa CKOPOCTH YacTHL *XUAKOCTH M0 rryOuHe BoOeMma. 

O6cystcdenue u 3akiio4enue. [penioKeHHbii MeTO ABIAeCTCA BeCbMa 93*PPeKTHBHBIM JIA PeWIeHHA 3aa4 oO 
MOABWOKHBIX Harpy3Kax, J[eHCTBYIOWIHX Ha Je AHO MOKPOB BOJOeMa KOHeYHO! TryOuHbI. OH MO3BOIAeT pellaTb 3aauH 
0 WeHCTBUM Harpy3KH, JBMKyWelica 10 eqAHOMY MOKpoBy M0 cIOX%KHOM TpaexTopun. IosyuenHble pe3syIbTaTbl MOryT 
ObITb HCHOJb3OBaHbI JIA pacueTa HallpMKeHHA HM TepeMeLeHHH JIEMOBOTO MOKpoOBa Mp MpOKaAKe JIEOBbIX JOpor WH 
CTPOHTEJIBCTBE aIPOAPOMOB Ha JIBALY. 


K.ouespie cs10Ba: OecKOHeUHBII Te qHOU TIOKPOB, ABMWKylWaAcd Harpy3Ka, MIpOW3BOJIbHad TpaecKTOpHA, MepeMeHHadA CKOPOCTb 


BuaarogapHocru. ABTop BbIpaxKaeT OaroapHOCTB peueH3eHTaM 3a YKa3aHHble 3aMe4aHHA, KOTOPbIe NOSBOJIMJIN 
TIOBbICHTb KAYCCTBO CTATbH. 


Aisa wutTupospannsa. Tanadypqun A.B. PaspadoTka MeToga pellieHua 3ayqa4uu ePpopMaluu seqaHoro NOKpoBa Toy 
TeVicTBHeM MpOv3BONbHO ABMKylWelica Harpy3ku. Advanced Engineering Research (Rostov-on-Don). 2024;24(2):170-177. 
https://doi.org/10.23947/2687-1653-2024-24-2-170-177 


Introduction. The development of the polar areas of the World Ocean and a large number of freezing inland 
reservoirs drive the need to study the fields of displacement and stresses of the ice cover caused by the action of various 
types of loads. Numerous papers by domestic and foreign scientists are devoted to solving these problems. Previously, 
it has been found that the mechanical properties of ice depend on its temperature and water salinity. Much attention 
was paid to the development of numerical ice models that accurately reflected the interaction of ice and ideal 
incompressible fluid. In [1, 2], the smoothed particle hydrodynamics was used for this purpose, in [3, 4] — the method 
of discrete elements. In [5], ice was modeled by an elastic plate lying on the surface of a stratified fluid. Models that 
allow cracks were considered in [6, 7]. Models of ice strengthened with reinforcing elements were presented in [8, 9]. 
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At this, in some papers, an ice cover is considered as an elastic plate lying on the surface of a reservoir [10, 11]. At 
the same time, in [12], on the basis of the conducted research, it is concluded that in some cases, the properties of ice are 
best described by the Kelvin-Voigt rheological model with one parameter (damping time). Therefore, numerous 
researchers use a viscoelastic plate when modeling the ice cover [13]. In [14], nonlinear models were used to describe the 
properties of ice. 

In some articles, the effect of a moving load on the ice cover was considered. In [15], the impact of a mobile load on 
the ice cover in a frozen channel was studied. In [16], the action of a load with an impulsive movement on the ice cover 
was described. Paper [17] is devoted to the study of the load moving along a frozen riverbed. Here, the rectilinear motion 
of the load was investigated [18]. However, in real conditions, it is often needed to deal with a load moving in a more 
complex way. Therefore, the objective of this work is to develop a method for solving problems about the action of a load 
moving over an ice cover along a complex trajectory. This will provide for a more accurate investigation of the effect of 
vehicles moving in a complex way on ice. 

This work is a continuation of research related to the problems about the effect of a moving load on various objects, 
whose results are presented in papers [19, 20]. 

Materials and Methods. Setting the task. A reservoir of finite depth with an infinite ice cover (an infinite plate), 
which is subject to the action of a vertical force moving in an arbitrary way — impulsively, is considered. It is assumed 
that the reservoir fluid is incompressible and executes irrotational motion. 

The problem is reduced to a system of differential equations [15]: 


x,y,t 
(14190, AW +07 O;W +kW +b0,F|_ = Olnrt) 
AF =0, 
where W(x,y,t) — ice cover deflection; E and 4 — Young's modulus and Poisson's ratio of ice, respectively; 


D = Eh3/12(1-7) — cylindrical bending stiffness; 4 — ice cover thickness; t. — strain relaxation time; A,” = (0,7+ 0,7)’; 
A= 0+ 0,?+0/; pxand ps— density of ice and water, respectively; c? = p,h/D; k = peg/D; b= ps/D; O(xy,t) — load 
acting on the ice surface; F(x,y,z,t) — speed potential. 

Under boundary conditions at z = 0 (ice-water boundary): 

OW =0.F. 
At the bottom of the reservoir at z = —H: 
0,F =0. 

In addition, it was assumed that the ice cover and the fluid in the reservoir were resting at the start time. The load was 
a concentrated unit force (one Newton) OQ (x, y, t), which moved arbitrarily along open free-form curve y. It was assumed 
that the displacement of the force was given in the form Q = Q(s(t)), where s — arc coordinate measured from some fixed 


= x(t) 


: x : 
point of trajectory y. The trajectory of movement was set parametrically in the form | ‘ where ¢ — time. 
= Volt 


The moving load was approximated by the expression: 
O(x,y,t) = 6° exp(-e" ((x —Xp (t))” +(y — Vo (1))” )\/e 


where  — numeric parameter. 
After applying the integral Fourier transform with respect to variables x and y, the integral Laplace transform with 
respect to ¢, formulas for calculating unknown functions W and F were obtained: 


W (x, y,t) mall 2p? /4s” Jo(pR(t-1))= ‘le ae O° — et" |dpdr, 


Pesan) asin [per tsa pete 


00 


to 
W (x,y,t)= xo) fren (pR(t -t)) 
00 


a [ete olata)s Japan, 


l(a +d)e (ar -(a - dye" |dpdt, 
R? (t)=87 +B?,5=x)(t) x,B=yo(t) J» 
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5 
ToP 


2(c*p + beth( pH)) : 


y= [ tp" -4(c p? + bpcth( pH))( p* +k] ,a= 


es ae 

2(c*p +bcth (pH)) 
Using the known relations from the theory of thin plates and the theory of the potential ideal flow, it is possible to 
obtain relations for calculating displacements and stresses in the ice cover, as well as the velocity vector component of 


fluid particles. 
When calculating the improper integral through numerical methods, the approximate relationship 


[, fle) did not 


io) A 
| f ( p) dp = | f ( p) dp was used, in which value A was chosen so large, that the error estimate 
0 0 


exceed the set value. 
Thus, for the amount of ice deflection 
tA € 


to ioe) 
1 1 1 
W (x,y,t) = wap | (1, p)dpdt = ral Us (+ P)drdre | fu (t, p)dpdt, 


this estimate has the form: 


t 0% 


1 2A? 2 
| Uo(t,p)dpdt Beg 
0A 


2nD nDy (A) 


y(A)=[ 13.4" -4(c?.4? + bActh( AH))( 44 +k] , 


Similar estimates can be obtained for other calculated quantities. These estimates were used to determine value A. 


2 
ef /40? 


—— did not exceed 
mDy(A) 


In the calculations carried out, value A was chosen such, that the estimate 


ten 

ssl Uo(t, p)dpdt| 0.001. 

When calculating the repeated integral, Simpson's quadrature formula (for variable t) and Chebyshev quadrature 
formula with equal weights for two nodes (for variable p) were used. The other values were calculated in the same way. 

Research Results. A method has been developed for solving problems on the action of a load moving along the ice 
cover of a reservoir filled with ideal fluid along a complex trajectory with variable speed. Using this method, calculations 
were carried out that showed the degree of impact of various parameters on the deformation of the ice cover. 

The described method does not impose restrictions on the shape of the trajectory of the concentrated force. In the 
calculations, a special case of a trajectory consisting of arcs of circles was considered (Fig. 1). The red dot indicates the 
position of the concentrated force at the time under consideration and the direction of movement of the force. 


Y, mo - - ~ ras “ = a 
10 | 
6% e 
2} 
| 
. i i i o i 1 4 n al 
0 4 8 12 16 X,m 


Fig. 1. Trajectory of concentrated force 
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The following parameter values were taken into account in the calculations: ice cover thickness h = 0.25 m, Young's 
modulus E = 500,000,000 H/m, Poisson's ratio of ice 1 = 1/3, ice density p = 900 kg/m’, fluid density p = 1,000 kg/m’, 
¢ =5. The calculation results are presented below. 

Figure 2 shows the change in the deflection of the ice cover at speed of force v = 2.5 m/s, tangential acceleration 


w; = | m/s”, reservoir depth H = 25 m, and relaxation time t) = 1s. 
The law of force motion along the trajectory was taken as: 
§=a,t? + aot? + ay. 
Coefficients a1, a2, a3 were selected in such a way that the force, being at the same point of the trajectory, had the 


required speed and tangential acceleration. 


W,m-107 - 


-15 4 


100 


Fig. 2. Change in ice cover deflection 


At other values of these parameters, the qualitative nature of the distribution of ice cover deflections remained 
almost unchanged. 
Z,m 


Y,m 


Fig. 3. Fluid motion caused by action of moving load on the ice cover 
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The fluid motion caused by the action of a moving load at the same values of velocity, tangential acceleration of the 
load movement, relaxation time, and depth of the reservoir is shown in Figure 3 (the distribution of the velocity vector of 


fluid particles is shown). 
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b) 
Fig. 4. Change in amount of maximum deflection of ice cover depending on: 


a — tangential acceleration rate; b — depth of reservoir 


The effect of the tangential acceleration of the force movement on the maximum deflection of the ice cover is shown 
in Figure 4 a. In this case, the force speed was equal to v = 17.5 m/s, and the relaxation time 1) = | s. 

Figure 4 b shows a graph of the dependence of the maximum deflection of the ice cover W on the reservoir depth H. 
Here, the speed of the load movement was v = 17.5 m/s, tangential acceleration w; = 1 m/s”, and relaxation time T = Is. 


W, m-10% 
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—1.4 


0 10 20 30 V, m/s 


Fig. 5. Change in maximum deflection of ice cover 
depending on force speed 


The dependence of the maximum deflection of the ice cover on the force speed is graphically shown in Figure 5. The 
depth of the reservoir was assumed to be 25 m, and the tangential acceleration — w,; = 1 m/s’. The solid line shows the 
dependence corresponding to the relaxation time To = | s, the dotted line corresponds to the relaxation time T = 10s. 
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Discussion and Conclusion. The influence of the reservoir depth on the maximum deflection of the ice was studied. 
A picture of the deflection of the ice cover was obtained due to the action of a concentrated force moving along a complex 
trajectory with variable speed. Calculations showed that with increasing depth of the reservoir, the maximum deflection 
of the ice cover decreased (Fig. 2). At the same time, a noticeable dependence of the deflection of the ice cover on the 
depth of the reservoir H occurred only for H<25 m. At great depths, the amount of the maximum deflections stabilized 
near a certain constant value and practically did not change. Thus, if H>25 m, then the depth of reservoirs can be 
considered infinite when calculating. 

An increase in tangential acceleration caused an increase in the deflection of the ice cover. Moreover, the dependence 
of deflection on tangential acceleration was very close to the linear dependence (Fig. 4). 

At low relaxation times tT, the speed of the load movement affected significantly the amount of ice deflection. At 
large times, the impact of the load movement speed on the deflection of the ice cover was noticeably reduced (Fig. 5). 

To study the state of the reservoir fluid, the distribution of the velocity vector of the fluid particles due to the action 
of the moving force on the ice was determined (Fig. 3). 

The developed method of solving problems and the results obtained with its help can be used in the construction of 
ice roads, design and construction of runways on ice. 
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